Electron Tunneling in Monolayer and Bilayer Graphene 
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I. INTRODUCTION 



Di Wu, Weiqiang Chen and Fu-Chun Zhang 
Department of Physics, and Center for Theoretical and Computational Physics, 
The University of Hong Kong, Pokfulam Road, Hong Kong, China 
(Dated: August 22, 2008) 

Electron's tunneling through potential barrier in monolayer and bilayer graphene lattices is in- 
vestigated by using full tight-binding model. Emphasis is placed on the resonance tunneling feature 
and inter-valley scattering probability. It is shown that normal incidence transmission probabilities 
| for monolayer and bilayer graphene exhibit different properties. Our calculation indicates that val- 

leytronics in graphene systems may be detected, generated and controlled by changing the structure 
parameters of the external electric potential. 

as: 
< 

Recently, great interest has been aroused in research on the physical properties of graphene, a onc-atom-thick al- 
lotrope of carbon, due to its successful fabrication in 2004—. Monolayer graphene is a truly two-dimensional material, 
with unusual electronic excitations described in terms of Dirac fermions that move in a curved spaced. The electrons 
in graphene seem to be almost insensitive to disorder and electron-electron interactions and have very long mean free 
patbA Hence, graphcnc's transport properties are rather different from what is found in usual metals and semiconduc- 
tors. Interestingly enough, these properties can be easily modified with the application of electric and magnetic fields, 
the addition of layers, and by controlling its geometry and chemical doping^. Apart from the interesting fundamental 
physics in this new system, graphene is attracting attention as a promising new material for electronic applications. 
For a review concerning the history, fabrication, fundamental properties, and future applications of graphene, we refer 
to the recent article^. 

The low energy charge carriers in graphene are described by a massless Dirac equation and have a linear energy 
dispersion which is isotropic near the Dirac points where the valence and conduction bands meet each other—. Such 
characteristics offer exciting opportunities for the occurrence of new tunneling phenomena and the development of high 
quality devices. Therefore, it may be valuable to investigate the electronic transport properties of graphene. In this re- 
gards, much of the phenomena associated with tunneling in graphene systems has been theoretically studieo— ^ i 10 ' 11 . 
O . It is interesting that, owing to the chiral nature of the quasiparticles, the propagation of charge carriers in monolayer 
graphene mimics the tunneling of massless fermions. This relativistic effect provides us an experiment test for the 
_— , Klein parado x 7 i 12 , which predicts that electron can pass through a high potential barrier to approach the perfect 
^ ' transmission. In contrast, for conventional non-relativistic particles, the transmission probability exponentially de- 
\ cays with the increasing of the barrier height. Besides this relativistic transport feature, other promising tunneling 
CO ■ properties of graphene systems are the ability to tune the carrier density through a gate voltage^, the absence of back 
scattering^, and the fact that graphene exhibits both spin and valley degrees of freedom^, which might be harnessed 
. . in envisaged spintroni o 15 ! 16 or valleytronic devices^. 

However, in these previous studies, considering the carriers as massless fermions remains an approximation, and we 
expect a deviation from a linear dispersion for high energies of the Dirac conei^. To this end, it is both important 
and interesting to study the transport properties in graphene system without this approximation. In this work, we 
investigate the tunneling properties of both monolayer and bilayer graphene lattice by using the full tight-binding 
model. This paper is organized as follows. In Sec. [Til we introduce the model of the system that we are considering 
and give explicit expressions for the wave functions in different regions through tunneling process. We also provide 
some details on how to compute the transmission amplitudes. In Sec. IIII1 we present the results for transmission 
probabilities in different barrier setups and the discussions. A brief summary and conclusion of the paper can be 
found in Sec. HVl 



II. MODEL 



We consider two kinds of lattice structure, each consisting respectively of monolayer graphene or bilayer graphene. 
A schematic picture of the system is shown in Fig. [TJ The basic structure involves a graphene sheet and a one 
dimensional trapezoid shape potential V(r) = V(x), which is y— independent. In the following we will introduce the 
tight-binding model of graphene with this external electric potential. 
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FIG. 1: Schematic of the external potential profile in graphene lattice. The energy of incident electron is Ei and the potential 
height is Vq. There are five different regions in this tunneling process which are marked by roman numbers. 



A. Tunneling in monolayer graphene lattice 




FIG. 2: Schematics of the lattice structure and the energy spectrum of monolayer graphene. Top left: Lattice structure of 
monolayer graphene, made out of two interpenetrating triangular lattice(ai and a2 are the lattice unit vectors, and Si, i = 1,2,3 
are the nearest neighbor vectors). Top right: The solid hexagon is a Brillouin zone. Dirac cones sit at the K+ and K- points. 
The reciprocal lattice vectors are bi and b-2. Bottom: The energy spectrum of electrons at a finite k y . There are four degenerate 
states for a given k y when incident energy Ei is small. Two of these states(ki and ki) have pseudo-spin +1, the other two 
states(k2 and ka) have pseudo-spin — 1. The solid arrow is the intra-valley scattering process(flip pseudo-spin) while the dashed 
arrow is the inter- valley scattering process(does not flip pseudo-spin). 

The honeycomb lattice of monolayer graphene can be described in terms of two triangular sublattices, A and B(see 
Fig. [2]). A unit cell contains two atoms, one of type A and one of type B. The lattice vectors can be written as 
ai = a(l/2, V3/2),a2 = a(— 1/2, y/3/2), where a sa 2.46A is the lattice constant. The reciprocal lattice vectors are 
given by b : = 2ZL(1, l/y/5) and b 2 = ^(1, -1/^3). 

In monolayer graphene, an atom of type A is connected to its nearest neighbors on B sites by three vectors Si. The 
nearest hopping tight-binding Hamiltonian describing this system has the forrn^ 



H = -t J2 (4 + h-c) , 



(1) 
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where i(« 2.8eV) is the nearest neighbor hopping energy, a a ^ and b a .j are the annihilation operators of electrons 
with spin a(a =Tji) on A and B sublattices, respectively. In momentum representation, the Hamiltonian reads 
#o = Ek, CT V4 ( k ) fto^V (k), where 



0(k) 

with (k) = — t e -' k(5 i = — e (k) e «¥>(k)_ This Hamiltonian acts on a spinor representing the wave function on two 

sublattices: ip a = (i/va, Vvb) T - Accordingly, quasiparticle's excitation spectrum has two branches(bands) with the 
dispersion E ± (k) = ±e (k) and^ 



fc^a yfikyCi k x a , . 

e(k) = ty 1 +4cos cos — ^-+4cos 2 -^-. (3) 

The upper band(i? + ) and the lower band(-E - ) meet at six corners of the first Brillouin zone. Only two out of the six 
are independent, which we choose to be K± = ±—(2/3,0) as shown in Fig. [51 With an expansion around K±, one 
can find that the Hamiltonian (2) reduces to a massless Dirac Hamiltonian with linear dispersion E^ (k) = ±frup|k|. 
It is important to note that in monolayer graphene the pseudo-spin is defined as eigenvalue of the Pauli matrix o\ . 

In tunneling problem, we consider the barrier of which the geometry is shown in Fig. [TJ The dynamics is governed 
by the Schrodinger equation for an incident particle of energy E, 

[H + V{x)]^ = E^. (4) 

We solve eqn. (j4]) by using the standard method. At first, we assume the incident electron wave propagates along 
the x axis with a given energy Ei and wave vector fcj. Because V has no y dependence, the wave vector along y 
direction is conserved, so the wave vector along the y-direction in all the regions remains to be ki y . In region I and 
V, the eigenstates with e(k) = Ei and k y = ki y are fourfold degenerate as shown in Fig. [2l while the corresponding 
eigenfunctions are 

^(kO = ^Ljw P V kir - (5) 



V2 \<t>W/Ei 

Inside the barrier, i.e, in region III, the electron's wave vectors should satisfy e(q;) = Vo — Ei,l = 1,2,3,4 and 
Qiy = kiy. The eigenfunctions inside the barrier can be written as 

After constructing all these wave functions, the next step is to determine incident and reflected waves through the 
whole tunneling process which can be done with the velocities of these states 



v x (k) 



at 



1 de (k) 

h dk x 

sin ( V) 



, \/3k a \ ( ka 

COS 7T— + 2 COS -S- 



2 




(7) 



Therefore, for the right movers in Fig. [TJ we require 2 ^ v x (k) > 0; on the contrary, for the left movers, their velocity 
v x (k) < 0. 

Before proceeding to the numerical calculation, it will be instructive to make some general analyses about the two 
different scattering processes in this system as shown in Fig. [2l In that figure, the solid arrow represents intra-vallcy 
scattering, while the dashed arrow represents inter-valley scattering. The latter one has a much larger momentum 
transfer (~ 2|K + |), hence it is usually neglected for low barrier (Vo <C t) tunneling ^ 10 ' 20 ' 21 . Furthermore, these 
two scattering processes are associated with the two operations on the pseudo-spin of the quasiparticles: intra-valley 
scattering flips the pseudo-spin; while during inter-valley scattering the pseudo-spin is conserved. 
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We can write down the general solutions in different region s 7 ' 10 . In region I, we have one incident wave (ki or 143) 
and two reflective waves (k 2 and k^), and the solutions can be expressed as 

*i(r) =?A(k 1 )+r 1 V'(k 2 )+r 2 ^(k4). (8) 

Here r\ is the intra-valley reflection amplitude and r 2 is the inter-valley reflection amplitude. In region III we have 
four hole states with 

*ni (r) = M (qi) + M (<fc) + M (<te) + M (q*) . (9) 
Finally, in region V, we have two transmitted waves 

*v(r) = M(ki) + *3^(ks). (10) 

Now we need to solve the Schrddinger equation in region II and IV numerically to find the wavefunctions. Upon 
applying the continuity of the wave function at the boundaries, one may obtain the coefficients r, /, and t. After 
that one can determine the transport probability by T = \t\\ 2 + \t 2 \ 2 X v x (k 3 ) jv x (ki). An important advantage of 
our model is that we can choose the incoming electron from different valleys (ki or k.3), which would lead to different 
tunneling properties in the results. 

B. Tunneling in bilayer graphene lattice 




FIG. 3: Schematics of the lattice structure and the energy spectrum of bilayer graphene. Top left: Interlayer coupling t± forms 
'dimers' from pairs of A2-B1 orbitals(dashed lines), which leading to the formation of high energy bands. Top right: Top view 
of bilayer lattice. The only hopping mechanism between two layers is Al(open)^B2(solid) hopping via A2-B1 (hashed) dimer 
state. Bottom: Band structure of bilayer graphene. At the Dirac point, the spectrum has a gap of t±. The pseudo-spin of 
lower band is +1, of upper band is —1. 

Bilayer graphene consists of two monolayer stacked as in natural graphite(see Fig. [3]). This so-called Bcrnal stacking 
yields a unit cell of four atoms(Al, Bl, A2 and B2) resulting in four electronic bands. The tight-binding Hamiltonian 
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for bilayer graphene can be written a a 5 ' 22 

H =-t ( a t,mi b <r,mj + h-C.) 

(i,i),m,tr 

-tx Yl (4,lA,2j + h.c.J , 
(U,2j),<r 



(11) 



where ij_(~ 0.35cV) is interlayer hopping through dimer states, m = 1,2 is plane index. In momentum space, the 
effective bilayer Hamiltonian has the form of 



Ha = 



( 





0* (k) 
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0(k) 






(12) 



The eigenstates of Eq. (|T2"j) are four component spinors ip a = (ip a Ai> V'crBi; VVA2) Vvb2) T j where V , ai,bi(V'A2,B2) are 
the envelop functions associated with the probability amplitudes at the respective sublattice sites of the lower(upper) 
graphene sheet. 

The Hamiltonian TCq determines the following spectrum of electrons in a bilayer graphene. There are four valley- 
degenerate bands, Eg, s = ±1, with 2 ^ 



£f(k)=± 



l + e 2 (k)+ S ^ 



(13) 



where e (k) is the dispersion of a monolayer graphene. The dispersion E^ t describe low energy bands while E^ 
describe higher energy bands as shown in Fig. [3j One can see that low energy excitations exhibit parabolic dispersion, 
while for larger k values, the linear E — k behavior is recovered. 

An important difference in the eigenfunctions between the monolayer and the bilayer graphene is that in the latter 
case there are eight eigenstates (two bands) for a given energy Ei and fixed k y . Accordingly, for four component 
spinor ip a , the pseudo-spin in bilayer graphene is defined as eigenvalue of Dirac matrix 



7i 



/0 1> 

10 

10 

Vl 0y 



(14) 



As a result, the pseudo-spin of state Ef in bilayer graphene is qps. So the lower band (E^) in Fig. [3] is associated 
with pseudo-spin +1, and the upper band (-Ei) with —1. Therefore, both inter- and intra- valley scattering do not flip 
the pseudo-spin. This feature is totally different from monolayer structure. Moreover, in our computation, the energy 
of incident wave satisfies JSj -C t± as this will likely be the experimental situation. So in region I and V we have four 
of these states(ki i 2,3,4) correspond to propagating waves(i?j^ 1 ) and the other four(/«i.2,3.4) to evanescent ones(-Ei). 
These evanescent modes have an complex value of the momentum in the a;— direction and must be considered to fulfill 
the boundary conditions. It is worth noting that the incident electron in bilayer lattice could only reside on band E'^ l 
with pseudo-spin +1. 

By solving Hamiltonian in Eq. (TT2")) . one finds the eigenvectors^ is external potential), 



$f (k) = C 



/a (k) <f>* 

( <* 



(k)/(^±- 
a(k) 

-V)/0(k) 
1 



(15) 



where C is normalization constant, and a (k) is defined as 

a(k) = -\(Et-V) 2 -\ct>(k)\ 2 ] /^(k). 



(16) 



Note again that for propagating waves, k x is real; while for evanescent solutions, k x is complex. 

Similar to the case of the monolayer graphene, by calculating the group velocity v x (k) = de (k) /hdk x , we can prove 
that ki.3 arc right movers (transmission waves) while k2.4 are left movers (reflective waves). On the other hand, to 
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select appropriate evanescent states, since Ki x (l — 1,2,3,4) are complex numbers, we should consider the asymptotic 
behavior of these states at ±00. 

Finally, the general solutions for bilaycr graphene can be expressed as^(incident electron coming from ki) 

*i (r) = (ki) + (k 2 ) + r 2 ^ a (k 4 ) 

+ r 3 V>+i (K2) + r 4 ^+i (ka) , (17a) 
4 

*hi (r) = ]T [frf-i + 9i*+i M . (17b) 

*v (r)=tiVii (ki)+t a ^! (ka) 

+ ta^+i + (k 3 ) • (17c) 

Here q; and r; are the corresponding wavevectors for propagating states (EZi) and evanescent states (-E+i) inside 
the barrier. Then by solving Eq. (j4|) numerically in II and IV, and utilizing continuity of the wave functions on four 
sublatticcs at boundaries, one may obtain the transmission coefficients for a bilaycr graphene lattice by using the 
same equations as in the case of monolayer structure. 



III. NUMERICAL RESULTS AND DISCUSSION 




FIG. 4: k y dependence of transmission coefficients in monolayer graphene for abrupt potential barrier(d = 0). The incident 
electron comes from valley K+ with energy 80meV. The potential height is 200meV with different widths D = 150a or 300a. 
During the whole tunneling process, we always have T + Rintra + Rinter — 1 to ensure current conservation. Note that different 
scaling are used in Rintra and Rinter 

We first calculate the transmission probabilities of charge carriers through monolayer graphene lattice. The results 
are depicted in Figs. [4] and [5] Fig. [4] shows examples of k y dependence of transmission probability for an abrupt 
potential barrier with height Vq = 200meV. Under this potential barrier, incident electron continues propagating as a 
hole in region III. The solid lines and dashed lines correspond to the potential width D = 150a and 300a, respectively. 
It is seen from the figure that, averagely speaking, the intra- valley scattering coefficient is much larger than inter- valley 
scattering coefficient. This is because in scattering process(see Fig. [2]), the intra-valley momentum transfer is much 
smaller than the inter-valley momentum transfer^ 2 |K + |), thus should have larger possibilities. It is also clear from 
the figure that the barrier remains nearly perfectly transparenfc£ii£(T ~ 1) for small k y . This phenomenon is unique 
to relativistic quasiparticles that incident electrons can be scattered into hole states inside the barrier. In this figure, 
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the incident particle is in state ki which comes from valley K + . We can also choose the incident particle coming 
from valley k 3 . The difference of these two results is very small, which means that valley discrepancy is very tiny in 
low-energy case. 
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FIG. 5: Transmission probability T and inter-valley scattering amplitude Rinter for normally incident electrons through mono- 
layer graphene as a function of external potential height. The incident electron comes from valley K+ with energy 80meV. The 
width of the barrier is d — and D = 150a or 300a. For normal incident tunneling, we always have Rintra = 0. 

Now let's examine the normally incident cases, i.e, k y = 0. The transmission coefficients are plotted in Fig. 
where we find that intra- valley scattering is totally suppressed. This can be understood in terms of the conservation 
of pseudo-spin. In normal incident tunneling, pseudo-spin becomes a good quantum number through the whole 
process. Therefore, spin flipping is not allowed when quasiparticlcs propagating in this system. Hence, inter-valley 
scattering is the only reflection mechanism. However, since momentum transfer kx — > k 4 is quite large, Ri n t er is very 
small (~ 0.01). So if we neglect the inter- valley scattering probability^, then monolayer graphene can be regarded 
as a condensed matter version Klein parado x 7 i 10 . Furthermore, we can see that the transmission amplitude shows a 
resonant feature that the envelope of T deceases monotonically as Vq increases. This can be understood as Dirac- 
like particle's resonance tunneling, where the resonance frequency and magnitudes varies significantly with different 
barrier widths. With increasing D, the resonance frequency increases obviously as shown in Fig.[5J while the envelope 
of transmission amplitude remains the same. On the other hand, the magnitude of oscillation also depends sensitively 
on barrier width d. By increasing d, inter-valley scattering amplitude falls noticeably. For example, in case of d = 3a, 
the inter-valley scattering amplitude falls to the order of 10~ 5 when Vb = 800meV. Thus one can achieve nearly 
perfect transmission even with very high barrier height by smoothing the potential step. 

By now we have investigated low-barrier (Vb <C t) tunneling in monolayer graphene by choosing incident particle 
as ki which comes from K + . If the incident electron comes from valley K_ , our result shows no obvious difference 
with valley K + . To amplify this valley-contrasting feature, we have to look at the high-barrier (Vb ~ t ~ 2.8cV) limit 
of the tunneling problem. The results of normally incident transmission arc depicted in Fig. [5] 

In this figure, one can see that the situation is completely different from low-barrier scattering. Besides the resonant 
feature in the transmission pattern, incident electrons coming from K + and K_ have totally distinct tunneling 
properties. The origin of this valley dissimilarity is due to the unique energy dispersion in graphene as shown in 
Fig. [21 At first, it is worth noting that pseudo-spin is a good quantum number in normally incident tunneling and the 
incident electrons might have two types of pseudo-spin, +1(K+) or — 1(K_). For low potential barrier, i.e. Vb < Ei+t, 
there are four propagating states in region III. Two of these states have pseudo-spin +1, while the other two have 
pseudo-spin — 1. Therefore, no matter which incident state we select, the pseudo-spin inside the barrier can always 
match the incoming quasiparticlc. So in low-barrier limit, resonant pattern is found in the transmission amplitude of 
both valleys. 

On the other hand, for high potential barrier, i.e. Vb > Ei + t, two propagating states with pseudo-spin —1 in 
region III turn to evanescent states^, while the propagating states with pseudo-spin +1 remain the same. As a 
result, different incident electrons come from different Dirac points K + and K_ must have distinct transmission 
properties. For incident electron in state ki with pseudo-spin +1, the potential barrier is always transparent. On 
the contrary, for electron coming from state k 3 with pseudo-spin —1, it cannot penetrate the potential barrier if 
Vb > Ei +t. Moreover, it is also shown in Fig.[S]that the possibility of inter-valley scattering can be notably enhanced 
if we increase the potential height Vb- Using these tunneling properties, one may obtain valley polarized electron 
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FIG. 6: Transmission probability T and inter-valley scattering amplitude Rinter for normally incident electrons through mono- 
layer graphene as a function of external potential height. The incident electron is of energy 80meV and may have valley 
polarization of Dirac points K+ or K_. The potential width is d = and D — 15a. Top: Incident electron comes from K+. 
Bottom: Incident electron comes from K_. 



current a 14 ' 17 in monolayer graphene by adding appropriate external gate voltage. 

Similar to low-barrier limit, we find that the magnitude of inter- valley scattering amplitude also depends very 
sensitively on potential height Vq and the thickness of region II, IV in Fig. [1] If we smooth the potential step by 
increasing d, then the envelope of inter-valley scattering coefficient shrinks much smaller very quickly. Furthermore, 
although inter-valley scattering requires a very large momentum transfer, we can still achieve very high scattering 
probability by increasing the barrier height. Our result shows that charge carriers in monolayer graphene, when 
subjected to a high potential barrier, exhibit totally distinct tunneling features from Dirac fermions. 

For the case of bilayer graphene, the computational results arc shown in the following figures. Some representatives 
k y dependence are depicted in Fig. \7\ and these should be contrasted to the case of monolayer in Fig. 5] The result 
is totally different since the bilayer energy spectrum(see Fig. [3]) splits into two bands. One can see that for potential 
barrier (Vq) higher than the energy of incident quasiparticle (Ei), it remains nearly totally reflective for small k y . In 
normally incident (k y = 0) case, it becomes totally reflective even though there are plenty of electronic states inside 
the barrieri2£. This behavior, which is similar to the tunneling property of non-chiral massive quasiparticles, is in 
obvious contrast to monolayer graphene, where massless Dirac fermions are always perfectly transmitted for small 

To understand this feature, let's take a look at the result of normally incident case as shown in Fig. [8] One can see 
in this figure that the transmission probability T ~ if 

E t <V <E l +t ± , (18) 

and the transmission amplitudes exhibit nonrcsonant feature within this range. This result is valley independent and 
can be explained in terms of pseudo-spin conservation in normal incident tunneling. As one can see in Fig. [3l under 
the condition of Eq. (|18[) . all the propagating states in region III are hole states EZ\ with pseudo-spin — 1. While 
the incoming wave, no matter which valley it comes from, must have pseudo-spin +1. As a result, the only survival 
tunneling mechanism is via evanescent hole states ET^. Therefore, (Ei,Ei becomes a non-penetrative region 

for bilayer graphene. This conclusion is coincide with the k y dependence as we depicted in Fig. [7l 

On the other hand, beyond this non-penetrative range, the resonant feature is found again for all barrier heights 
as the pseudo-spin of propagating states in region III can always match the incoming electron. Further more, valley 
difference is not obvious in low energy limit when V <C t, which is the same as the case in monolayer graphene 
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FIG. 7: k y dependence of transmission coefficients in bilayer graphene for abrupt potential barrier(d = 0). The incident electron 
comes from valley K+ with energy 80meV. The potential height is 120meV with different widths D = 150a or 300a. 



under low-barrier limit. While for high barrier tunneling, the transmission patterns of different Dirac points are 
notably distinct. This distinction is due to the inter- valley momentum transfer |ki — ki4 1 = 2|ki| is larger than 
|k3 — lc2 1 = 2 1 lea | So intrinsically the probability of inter- valley scattering amplitudes ki — > k4 should smaller than 
k,3 — > k2. In Fig. [8l this difference is amplified by high potential barrier and shows the existence of valley polarized 
current in bilayer graphene^. Moreover, in Fig. [5J the patterns of reflection coefficient Rintra are almost the same. 
This can be understood as intra-valley momentum transfer in two valleys are equal to each other 

|kx -ka| = |k 3 -k 4 |. (19) 

Finally, we also point out that the dependence of transmission amplitude on d and D in bilayer graphene is the 
same as the result in monolayer structure. If we increase the width of region III, the envelope of T curve does not 
change while the period of resonant oscillation decreases. On the other hand, if we smooth the potential step by 
broadening the width of region II and IV, the amplitude of intcr-vallcy scattering falls very quickly to smaller than 
10 -5 . However, intra-valley scattering probability does not depend on d in our computation. Thus we can tune the 
proportion of these two scattering mechanisms separately by adding arbitrary electronic fields at the edges of the 
external potential barrier. 

Before closing, just a comment on our calculation. Actually, previous studies of tunneling properties in monolayer 
(bilayer) graphene that were based on a continuum model have used a massless Dirac fermion (massive excitations) 
approximation, which is expected to be accurate for low-energy quasiparticle3^'^ i 10 ' 20 i 22 . Hence, inter-valley scatter- 
ing is rarely considered in investigating the transport properties of graphene systems 2 ^ because of the large separation 
of Dirac points in momentum space. Moreover, valley contrasting physics is also absent in low-energy quasiparticlc's 
tunneling features in this system. However, our result shows that, although low-energy tunneling in monolayer (bi- 
layer) graphene is sufficiently described by massless (massive) fermion approximation, in high-energy limit there will 
be a significant deviation between this approximation and the real system^. With the increasing of barrier height, 
inter-valley scattering can be realized and amplified significantly as shown in Figs. [HI and [5] At the same time, we 
find that valley discrepancy in the tunneling problem would be gradually magnified by strengthening the external po- 
tential. These effects, which are distinct from the properties of low-energy limit, lead to the possibility of generating, 
detecting and controlling valley polarized currents in graphene systems by electric mean s 1 i 1 . 
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FIG. 8: Transmission coefficients for normally incident waves through bilayer graphene as a function of external potential 
height. The incident electron is of energy 80meV and may have valley polarization of Dirac points K+ or K_. Potential width 
is d — and D — 30a. Top: Incident electron comes from K+. Bottom: Incident electron comes from K-. 



IV. SUMMARY AND CONCLUSION 



Based on a tight-binding model, we have investigated single particle tunneling properties through potential barrier 
in monolayer and bilayer graphene lattices. The normal incidence and angular dependent transmission probabilities for 
two kinds of graphene structure have been numerically calculated. Our result illustrates that the tunneling behavior 
in graphene is much richer than what was anticipated in masslcss Dirac fcrmion approximation. Furthermore, it is 
shown that both of the intra- and inter-valley scattering probability are strongly dependent on the applied external 
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potential, forming the basis for the valley-based electronics applications of these systems. Our results may provide an 
important reference to the design of electron devices based on graphene materials. 
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